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INTRODUCTION 


MOST previous studies of hybridization have been restricted to its 
qualitative aspects. Quantitative measurements of the fertility rela- 
tionships within and between species are urgently needed. To meet 
this need some simple formulae are here presented which have been 
derived in an attempt to evaluate the fertility of parent stocks as well 
as the amount of reproductive isolation which is involved in intersub- 
specific and interspecific crosses in Peromyscus. 

Mice of the genus Peromyscus offer a unique opportunity for the 
study of hybridization. They occupy a variety of ecologic niches and 
are divided into many species and subspecies. Most species can be 
bred in the laboratory, and some interspecific and intersubspecific 
crosses are successful (Dice, 1933). 
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MATERIALS AND METHODS 


The stocks of mice used are as follows: 

Peromyscus maniculatus blandus.— Taken in the vicinity of Tularosa, 
Otero County, New Mexico, in 1938, by W. F. andF. A. Blair. 

P. maniculatus bairdi. — Taken in Washtenaw County, Michigan, in 
1946. 

P. polionotus leucocephalus. — Taken on Santa Rosa Island, Florida, 
in 1942, by W. F. Blair. 

The animals were reared by the general methods described by Dice 
(1947). Only individuals from about two months to one year of age were 
used in the matings. Close inbreeding was avoided, and all matings in 
which the relationship was closer than first cousin were excluded. Each 
mating was examined for litters every day for a period of six months, 
starting three weeks after the mating was made. Those matings which 
produced one or more live-born young were considered to be fertile. 
The sterility of some mated females was later confirmed by the ab- 
sence of placental scars in their uteri. 


Formulae for the Measurement of Fertility 
and of Reproductive Isolation 


Before an attempt is made to analyze or to interpret the experimen- 
tal results, certain formulae needed for the measurement of fertility 
and of reproductive isolation are presented and, when necessary, a 
brief explanation is given of their derivation. 

Only a fraction of laboratory-bred animals is fertile, as is probably 
true of field-living animals. This is undoubtedly due to many factors, 
involving both heredity and environment. It is first necessary, there- 
fore, to make an estimate of the probability of an individual's being fer- 
tile. If b is the number of fertile matings among the total number, n,, 
of matings made, and if it is assumed that f, the probability of fertility, 
is equal and independent in the two sexes (this point will be discussed 
later), then the proportion of fertile matings is equal to the probability 
of a fertile female being mated with a fertile male, or b/n, = f?. From 
this relationship the probability of an individual mouse being fertile is 
readily obtained as f =vb/n,. 

Therefore, the following relations are obtained: 

Proportion of fertile matings, 


pe | (1) 


ny 
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Probability of an individual being fertile, 


Ib 
fe =. (2) 


The standard deviation (S) of f is easily obtained by the maximum 
likelihood method (Mather, 1951), as follows: 


“of sl aif? (3) 
S, =/ ripe 


Since it is not infrequent that certain males have been previously 
proved to be fertile, the probability of these males being fertile is 
unity, and as a result the proportion of fertile matings, c, among total 
matings, n,, involving these males should be equal to the proportion of 
fertile females among the total number of females used. That is: 

Probability of an individual fenvale being fertile, 


c 
f-= =, (4 
Tl ) 
The standard deviation in this case is: 
{(1-f) 
S;= Pee (5) 


Consequently, two different f values, f, and f,, may be calculated 
from the separate sets of data obtained from different mating types, 
that is, f, is calculated from one set of data gathered from matings in 
which the fertility of both parents is unknown, and f, is estimated from 
another set of data gathered from matings in which the fertility of only 
the female parent is unknown. The significance of the difference of 
these two f values from zero can be tested by the standard deviation 
method. Thus, if the standard deviation of f, is S;,, and that of f, is 
S-,, and if the two sets of data are independent, then the standard de- 
viation of (f, + f,) is: 


Seice,= /97, + SG, 
where S,,and Sr, are shown in equations (3) and (5), respectively. 
Hence the difference between (f, -f,) and zero in terms of the standard 
deviation of (f, - f,) is: 


genie be 


¥ Si.+ SF, 
and the probability associated with any t value can be obtained from a 
table of the normal distribution. 
If these two f values differ statistically as tested by the standard de- 
viation method, then let f' be the probability of a female being fertile, 
and substitute it in equation (1): 


pee 
ny 
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or 
ee (6) 
in, 


Since b, n,, and f' are known, equation (6) is readily solved for f, 
which is the probability of an individual male being fertile. But f and 
f{' are not obtained independently and there are no means of testing 
whether the difference between f and f' is significantly different from 
zero. When these two f values are known, all the other equations with 
the exception of (7) and (8) can be readily solved. 

In my experiments, except for the matings of blandus, two separate 
sets of data are not available. It is assumed that f and f' are of the 
same magnitude, and therefore the prime of f will be dropped in the 
following equations. 

If there are two separate sets of data for any kind of matings, and 
if f can not be shown to differ in the two sets, it is logical to pool them 
to obtain a combined estimate of f by the maximum likelihood method. 

Now, the joint probability, P, of obtaining precisely b fertile mat- 
ings among n, matings when both parents are of unknown fertility and c 
fertile matings among n, matings when the male parent is known to be 
fertile is: 


b So = n, b oa irate Nn, c No-c 
P (7% )= [(pi) (Ee) (A -B)ms-4][(%) (Oe (1 - 1)", 
Intuitively, the best estimate of f is that value which maximizes P, 


that is, that estimate which makes most probable the observed sampling 
distribution. Now 


log P = log K, +blog f? + (n, -b) log (1-£?) + log K, + c log f + (n,-c) log (1-f) 


where 


and both are constants depending on the coefficients of the likelihood term. 
Maximizing this logarithm gives: 


d logub Se 2htcn(n abies. © are 
rYsuiwem aegis artis G 


Equation (7) is a quadratic equation and is easily solved for f. The 
standard deviation of the maximum likelihood statistic f (Mather, 1951) 


1S: 
te {(1-f?) s 
S:" /4n,.€ 4, et (8) 


When the probability of individual fertility of each parent stock is 
known, the reproductive isolation index of the intersubspecific and inter- 
specific crosses may be evaluated. Let f, and f, be the probabilities of 
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individual fertility for two parent stocks A and B, respectively, and 
f?,, the proportion of fertile matings among the total cross matings of 
the type female stock A times male stock B. If no disturbing factors 
are involved, such as psychologic or gametic isolation, one would ex- 
pect: 


Po sent feew 


The reproductive isolation index may now be defined as the ratio of 
the difference between the product of the two probabilities of individual 
fertility for two parent stocks and the observed proportion of fertile 
matings among the total cross matings to the sum of these two quanti- 
ties. Therefore, if one measures the reproductive isolation between 
the female stock A and male stock B, the formula is: 


re Rigs ES 
eel i ESTOS (9) 


Here f,f, is used in both the numerator and the denominator to ac- 
count for the effect of captivity in the laboratory. Although the labora- 
tory conditions may be made as suitable to the animals as possible, 
they are never perfect, and the fertility of one stock might be more 
influenced by captivity than that of the other. One would, therefore, 
expect the fertility of cross matings to be affected in varying degrees 
according to the parent stocks used. In other words, cross fertility 
is a function of the probabilities of individual fertility in the two parent 
stocks. If, therefore, the proportion of fertile matings in an outcross 
is equal to the product of the two probabilities of individual fertility of 
two parent stocks, the reproductive isolation index is equal to zero. 

In all other cases the values may fluctuate in either direction, but it 

is easily seen that the values are bounded between -1.0 and +1.0. If 

all cross matings are sterile, the value of the index is +1.0. If nega- 
tive values of the index are obtained, cross mating occurs in preference 
to matings between individuals of the same stocks. 

Similarly, to measure the reproductive isolation of the reciprocal 
cross, i.e., the mating between female stock B and male stock A, the 
equation is: 


tse ee Fe 
Ihe = FEE 10) 


Finally, these two types of reciprocal matings may be combined, 


and their joint reproductive isolation index, I,,,), measured The for- 
mula is: 
es fot, - fab) 11 
ToT ia Ey, (11) 
where f?_, is the proportion of fertile matings in both reciprocal crosses 


(ab) ; ; 
of stocks A xB. Inusing this formula, it is desirable that the number 


of matings in the two reciprocal crosses be equal. The effect on the 
reproductive isolation index of marked departures from equality in the 


6 Be ad eg 8 CE Vere 


number of matings in the two reciprocal crosses has not been explored. 

The statistical significance of departures of the observed values 
from the expected values, that is, f,f,, can readily be tested by chi- 
square as shown by Table I. 


TABLE I 


Scheme for Chi-square Test 


Fertile Matings Sterile Matings 


Observed b (b - n,f, f,)? 
n,f,f, (1 - f f,) 


Mating 


n, -b 


Cross between 
female stock 
A and male 
stock B 


ny (1 r f, f,) 


Expected n,f,f, 


Cross between | Observed c n-c (c - nf, £,)? 


err arr 1 
female stock f,f(1 - f,f,) 
B and male __| Expected mip n, (1 - f,£,) Mat 5 
stock A 
Both recipro- | Observed bic n,+n,-b-c (+c) - (n, +n,)£,f,)? 1 


calcrosses | Fenected (n,+n,)f,f, | (mn, +m,)(1- f,&) | (n,+n,)f,f, (1 - fal) 


Analysis of Chi-square 


Total Sum of two X?'s for two 2 
reciprocal crosses 
Deviation x’ from the item of both 1 


reciprocal crosses 


Heterogeneity Total X? minus deviation x 


From Table I two kinds of chi-square are obtained for comparison. 
One is the total chi-square, which is obtained from the sum of two chi- 
squares for two reciprocal crosses, and the other is the deviation 
chi-square, which is obtained directly from the item of both reciprocal 
crosses. The difference between these two chi-squares is the hetero- 
geneity chi-square, which is used to test the agreement of the two re- 
ciprocal crosses. The analysis is presented in tabular form. Details 
will be given later in the numerical analysis. 


ANALYSIS OF DATA 


The results of breeding within each of the three parent stocks, 


namely, P. p. leucocephalus, P. m. blandus, and P. m. bairdi, are 
summarized in Table II. 


In these parent stocks the probability of individual fertility is highest 
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in blandus. The probability value of blandus is 0.94, that of bairdi is 
0.69, and that of leucocephalus is low, only 0.56. The wide difference 
in the probability of individual fertility among these three parent stocks 
suggests that, among other things, the laboratory conditions were more 
favorable to the blandus stock than to the others. 

In the case of blandus, data from two separate sets of matings are 
available, namely, one set from matings in which the fertility of 
both sexes is unknown, and another from the matings in which the fer- 
tility of the males is known from previous matings. The fertility in 
these two sets of matings, as measured by their f values, is not sig- 
nificantly different statistically (t = 0.94, P = 0.40-0.30). The best 
estimation of the f value may be obtained by applying the maximum 
likelihood method (equation 7), with the following results: 

2(28) 2(4)f 18 3 


Mamie kee 
The solution of this equation gives an f value of 0.91, and substitution 
of this f value in equation (8) gives S; = 0.03. Since this probability of 
individual fertility and its standard deviation are obtained from the com- 
bined estimation of two sets of data, these values will be used for 
blandus in later calculations. 


TABLE III 


Reproductive Isolation Index of Certain Intersubspecific 
and Interspecific Crosses in Peromyscus 


Proportion of Reproductive 
Fertile Matings | Isolation Index 
Among Total (f?) (I) 


Sterile 


Mating : 
Matings 


@ P. p. leucocephalus 
x 
o P. m, blandus 


@ P. m. blandus 
x 
o& P. p. leucocephalus 


Both reciprocal 
crosses 


9 P. m, bairdi 
x 
o P. m. blandus 


9 P. m. blandus’ 
x 
o P. m. bairdi 


Both reciprocal 
crosses 
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Since the probability of individual fertility in the three parent stocks 
is known, equations (9), (10), and (11) are used to derive reproductive 
isolation indices for the intersubspecific and interspecific crosses. 
Their results are shown in Table III, and their chi-square values are 
shown in Tables IV and V, respectively. 


TABLE IV 


Chi-square Test for Reproductive Isolation in Crosses 
between P. p. leucocephalus and P. m. blandus 


-—— — 
Mating a Fertile Sterile Total - a df Pp 
Matings Matings x ee 
Se Ne — 
g FP. p. leucocephalus |Observed 10 54 64 
x SUNO TA) of <0.01 
o P.m. blandus Expected 32.61] 31.39 64 
¢ P.m. blandus Observed 10 21 31 
x 4.34) 1 0.05-0.02 
o& P. p. leucocephalus | Expected 15.80] 15.20 31 
Both reciprocal Observed 20 75 95 
crosses 33,99 | 1 <0.01 


46.59 95 


Expected 48.41 


Analysis of Chi-square 


Item | Ve | d. f. aot Pp 


Total 36.31 2 < 0.01 
Deviation 33.99 1 < 0.01 
Heterogeneity 2.32 1 0.20-0.10 


The reproductive isolation index between ?leucocephalus and &* 
blandus is quite high. Its value, 0.52, is statistically significant. The 
reproductive isolation index of the reciprocal cross is only 0.23, which 
is barely significant (P = 0.05-0.02). The joint reproductive isolation 
index of the two reciprocal crosses is 0.42, which, again, is statisti- 
cally significant. Finally, the heterogeneity test shows that the data 
from these two reciprocal crosses are homogeneous. 

The over-all picture given by these figures indicates that a certain 
amount of reproductive isolation has developed between leucocephalus 
and blandus, but that the degree of reproductive isolation is somewhat 
different in the reciprocal crosses. These two species occupy different 
geographic areas and ecologic niches with no naturally interbreeding in- 
termediates. Since there is no chance for gene exchange between them in 
nature, genetic differences presumably have accumulated, and gradually 
reproductive isolation has been built up. In addition, leucocephalus is 
smaller in body size. When this species is used as the female parent 
in the cross, difficulty frequently arises in parturition (Watson, 1942). 
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It has been found that over half of the leucocephalus mothers died with 
their fetuses, but such high mortality does not occur if the larger blandus 
is used as the female parent. This is one of the facts which give an ex- 
planation of the difference of reproductive isolation indices between 
these two reciprocal crosses. 

Reproductive isolation between blandus and bairdi presents differ - 
ent statistics. The reproductive isolation index of the cross between 
9 bairdi and o blandus is 0.35, which is statistically significant, and that 
of the reciprocal cross is -0.10, statistically nonsignificant, whereas 


TABLE V 


Chi-square Test for Reproductive Isolation in Crosses 
between P. m. blandus and P. m. bairdi 


— = 
Mati Fertile Sterile 
ee Matings Matings 
9 P. m. bairdi Observed 
x 21.91 <0.01 
o P. m. blandus Expected 29.51] 17.49 47 
@ P. m. blandus Observed 34 11 45 
x 3.13 0.10-0.05 
o P. m. bairdi Expected 28.26] 16.74 45 
Both reciprocal Observed 48 aa 92 
crosses 4.44 0.05-0.02 


| Expected 57.77] 34.23 92 


Analysis of Chi-square 


Item 


Total <0.01 
Deviation 0.05-0.02 
Heterogeneity < 0.01 


the joint reproductive isolation index of the two reciprocal crosses is» 
0.09, barely significant (P = 0.05-0.02). The chi-square test for het- 
erogeneity shows that the fertility is different in the two reciprocal 
crosses. 

These findings indicate that some reproductive isolation between 
bairdi and blandus probably exists, but that the joint reproductive iso- 
lation in the two reciprocal crosses is negligible (0.09). If one inspects 
the geographic distribution patterns of these two subspecies, it is found 
that intermediate interbreeding forms exist (Osgood, 1909). These 
facts coincide with the Laboratory findings that reproductive isolation 
between these two subspecies is negligible. But there is a significant 
heterogeneity, i.e., the fertility of the two reciprocal crosses is not 
the same. The question arises, therefore, of why a particular cross 
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between bairdi and blandus shows reproductive isolation but the recip- 
rocal does not. This is a difficult question to answer. Nevertheless, 
a clue may be found in the difference in individual fertility of the two 
parent stocks. It seems likely that the female may play a more im- 
portant role in reproductive isolation than the male. To demonstrate 
whether this is true or not will require further experiments. 


SUMMARY 


Several simple formulae are presented for the measurement of fer- 
tility within and between parent stocks. The probability of individual 
fertility (f) of a parent stock, on which all other calculations depend, 
is equal to the square root of the proportion of b fertile matings among 
the total n,, when the fertility of both sexes is unknown, or 


r=, 
ny 
and is equal to the proportion of c fertile matings among the total n,, 


when one sex is known to be fertile from a previous mating, or 


Cc 
n 


f= 
2 

The reproductive isolation index is used to measure the reproductive 
isolation quantitatively and is defined as the ratio of (f,f,-f%,,), the differ- 
ence between the product of the two probabilities of individual fertility 
for two parent stocks and the observed proportion of fertile matings 
among the total cross matings, to (f,f,+ f?,,), the sum of these two 
quantities, or ae fae, 

om Eta fg 2 


The statistical significance of the departures of the observed values 
from the expected values is ascertained by the chi-square test. 

The probabilities of individual fertility thus obtained are quite dif- 
ferent in the three parent stocks studied. It is highest in Peromyscus 
maniculatus blandus (0.91), lower in P. m. bairdi (0.69), and lowest 
in P. polionotus leucocephalus (0.56). This probably means that, 
among other things, the laboratory conditions were more favorable to 
blandus than to leucocephalus. 

The reproductive isolation indices of the cross leucocephalus x 
blandus and of its reciprocal cross are rather high (0.52 and 0. 23), 
hence their joint reproductive isolation index is also high (0.42) and 
statistically significant. This shows that some reproductive isolation 
exists between these two species. The reproductive isolation indices 
of the cross bairdi x blandus and its reciprocal are either low or nega- 
tive (0.35 and -0.10); in turn, the joint reproductive isolation index is 
very small (0.09) and statistically barely significant. This indicates 
that the reproductive isolation between these two subspecies is negli- 


gible. 
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